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ABSTRACT 

We have mapped the inner 360" regions of M51 in the 158 /mi [C II] line at 55" spatial resolution 
using the Far-infrared Imaging Fabry-Perot Interferometer (FIFI) on the Kuiper Airborne Observatory 
(KAO). The emission is peaked at the nucleus, but is detectable over the entire region mapped, which 
covers much of the optical disk of the galaxy. There are also two strong secondary peaks at ~ 43% to 
70% of the nuclear value located roughly 120" to the north-east, and south-west of the nucleus. These 
secondary peaks are at the same distance from the nucleus as the corotation radius of the density wave 
pattern. The density wave also terminates at this location, and the outlying spiral structure is attributed 
to material clumping due to the interaction between M51 and NGC 5195. This orbit crowding results in 
cloud-cloud collisions, stimulating star formation, that we see as enhanced [C II] line emission. The [C II] 
emission at the peaks originates mainly from photodissociation regions (PDRs) formed on the surfaces 
of molecular clouds that are exposed to OB starlight, so that these [C II] peaks trace star formation 
peaks in M51. The total mass of [C II] emitting photodissociated gas is ~ 2.6 x 10 8 M Q , or about 2% of 
the molecular gas as estimated from its CO(l— >0) line emission. At the peak [C II] positions, the PDR 
gas mass to total gas mass fraction is somewhat higher, 3— T7%, and at the secondary peaks the mass 
fraction of the [C II] emitting photodissociated gas can be as high as 72% of the molecular mass. Using 
PDR models, we estimate the far-UV field intensities are a few 100 times the local Galactic interstellar 
radiation field, similar to that found near OB star forming giant molecular clouds in the Milky Way. The 
density solution is degenerate, with both a low (n ~ 10 2 — 10 3 cm" 3 ), and a high (n ~ 10 3 — 10 6 cm" 3 ), 
density solution. Our analysis shows that a substantial amount of the observed [C II] emission from the 
galaxy as a whole can arise from the ionized medium, and that the contribution from the cold neutral 
medium (CNM) is not negligible. At the [C II] peaks, probably ~ 7 - 36% of the [C II] emission arises 
from the CNM, while north-west of the nucleus, most of the observed emission may arise from the CNM. 

Subject headings: galaxies: interaction, individual(M51) -infrared: galaxies 



1. INTRODUCTION 

We are carrying out a study of the properties of the 
interstellar medium (ISM) in galaxies to understand the 
effects of the interplay between the global conditions of a 
galaxy and its star formation activity. The global proper- 
ties of a galaxy most likely set the preconditions for star 
formation and also trigger the star formation. By "global" 
we are refering to the conditions in the spiral arms, in- 
terarm medium, nuclei, or individual regions in the few 
hundred parsec range. Interacting galaxies are excellent 
laboratories for this investigation because the interaction 
can speed up processes (e.g. increase mass flow rates) or 
enhance physical conditions (e.g. increase densities) in the 
interstellar medium. Since every galaxy is unique in some 
sense it is important to carry out case studies of individual 
interacting galaxies. Here we have chosen the galaxy M51 
to investigate its ISM and its star formation activity. 

The galaxy M51 (NGC 5194), also called the whirlpool 
galaxy, is a grand-design spiral of Hubble type Sbc. It has 
an inclination angle of w 20° (Tully 1974 b) and thus it is 
seen almost face on. The total spatial extent of M51, as 
seen in the visible, is about 7'. M51 is interacting with its 
companion NGC 5195 which is 4.5' to the north. Due to 



its proximity (9.6 Mpc, Sandage & Tammann 1975) and 
its face on appearance, M51 is one of the best studied in- 
teracting galaxies. Despite this fact some aspects of this 
galaxy remain a puzzle. 

The star formation rate of M51 is not spectacular. It 
is quite normal for an isolated Sbc-Sc galaxy (Kennicutt 
1998). However, M51 might have been a Sb galaxy prior 
to the interaction with NGC 5195 and evolved to a Sbc- 
Sc late-type galaxy (Tully 1974 c; Kennicutt 1998). The 
galaxy interaction is also believed to be responsible for 
the current density wave in M51 enhancing the inner spi- 
ral arm structure and triggering star formation in the spi- 
ral arms. The density wave pattern, however, terminates 
at corotation and the spiral arms beyond corotation are 
attributed to material clumping induced directly by the 
galaxy interaction (Tully 1974 c). It is clear from the ar- 
rangement of the ionized gas (Tully 1974 a,b,c; van der 
Hulst et al. 1988; Tilanus & Allen 1991), molecular gas 
(Garcia-Burillo, Guclin, & Cernicharo 1993; Lo et al. 1987; 
Rand & Kulkarni 1990; Vogel, Kulkarni, & Scoville 1988), 
atomic gas (Tilanus & Allen 1989, 1991; Rots et al. 1990), 
and dust in the spiral arms that star formation in the arms 
is due to a density wave. However, the crowding of H II re- 
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gions in the north-east and the south-west of M51, as seen 
for example in the Ha map of van der Hulst et al. (1988), 
may not be caused by the density wave but directly linked 
to the galaxy interaction. In this case the galaxy interac- 
tion has strong and direct control over the location and 
the strength of the star formation in M51. 

To investigate the star formation activity, we observed 
M51 in the 158 /im [C II] fine-structure line. The [C II] line 
emission from galaxies arises predominantly in the warm, 
dense, photodissociated surfaces of molecular clouds ex- 
posed to starlight from nearby OB stars (e.g. Stacey et al. 
1991). These regions, commonly called photodissociation 
regions (PDRs) are heated predominantly by photoelectric 
ejection of energetic electrons from grains, and cooled by 
far- infrared (FIR) fine structure line radiation from 0°, 
and C + . Typically, the strongest cooling line is the [C II] 
line, but the [O I] 63 fj,m line is often roughly the same 
strength in starburst nuclei. The [C II] line emission from 
moderate velocity shocks is much weaker than that from 
PDRs (Hollenbach and McKee 1989). 
2. OBSERVATION 

The observation of the [C II] 2 Pz/i^> 2 Pi/i nn e structure 
line at 157.7409 /im was carried out with the MPE/UCB 
Far-infrared Imaging Fabry-Perot Interferometer, FIFI 
(Poglitsch et al. 1991; Stacey et al. 1992), on board the 
Kuiper Airborne Observatory in June 1994 and April 1995. 
The field of view of the 5 x 5 array was 200" x 200" 
(40" x 40" per pixel) and the beam size was 55" (« 
8.3 x 10~ 8 sr; FWHM). We placed the array at 7 positions 
in M51. Four positions were placed 160" apart in a square 
to cover the whole galaxy. To fully sample the central re- 
gion, two positions were shifted diagonally in north-east 
and south-west directions by 1.5 pixels with respect to the 
center and one position was placed at the center. The 
center position ((0,0) position) is at RA = 13 h 27 m 46.4 s , 
Dec = 47°27'13" (1950). The total observed field of view 
was 360" x 360". The position accuracy is estimated to be 
about 10". 

Except for the position at the center all scans were car- 
ried out in "stare" mode (fixed plate separation of the 
resolution determining Fabry-Perot interferometer). The 
spectral resolution in this mode was 120 km s _1 (FWHM). 
To take the dispersion of the velocity shifts in M51 into 
account (more than 100 km s _1 ) we readjusted the plate 
separation of the Fabry-Perot interferometers for each in- 
dividual array position. However, the spectral range of 
some of the edge pixel which cover outer parts of M51 do 
not cover the entire expected velocity range of the [C II] 
line. Those pixels were not used for the data analyses. 
For the center position we scanned the Fabry-Perot inter- 
ferometer over 300 km s _1 . The spectral resolution at this 
position was 114 km s _1 (FWHM). We used the H 2 S line 
at 158.0148 jum to calibrate the transmission wavelength 
of the Fabry-Perot interferometers. 

Internal blackbodies were used to Hatfield the detector 
array. We observed the planets Jupiter (June 1994) and 
Mars (April 1995) for absolute calibration of our data. For 
the brightness temperatures of Jupiter at 158 /xm we used 
T = 128 K (Hildebrand et al. 1985). We determined the 
brightness temperature of Mars to be 207 K at 158 /im 
from the published values of Wright (1976) and Wright 
& Odenwald (1980). The mean diameter of Jupiter at the 



time of the observation was about 42" and the mean diam- 
eter of Mars was about 10". The accuracy of the absolute 
calibration is about 30%. 

The final map of the integrated [C II] line intensity (Fig- 
ure 1) was made using a maximum entropy image restora- 
tion program. This program was written especially for the 
FIFI data and follows the algorithm as outlined in Skilling 
& Bryan (1984). It takes the data from all positions ob- 
served with FIFI and reconstructs an image on a finer 
grid using the maximum entropy method and constrains 
the data points by a chi-square fit. The result is a decon- 
volved image which is then convolved with a 55" (FWHM) 
Gaussian beam. An additional large error term was added 
to pixels which do not cover the full velocity range within 
their beam. This resulted in a smoothing effect due to a 
strong influence from neighboring pixels which do cover 
the full vcloctiy range. The contour lines in the final map 
are smoother then in the raw map and all features can 
be traced back to the observed data (no artificial features 
were introduced by the program). On the other hand all 
obvious features in the observed data appear in the final 
map. 



3. RESULTS 

3.1. Morphology 

The contour map of the integrated [C II] line inten- 
sity superimposed on an optical image of M51 (Figure 1) 
shows that the [C II] emission is distributed over the 
whole galaxy, covering a total solid angle of 2.85 x 10~ 6 sr 
and peaking at the nucleus. The integrated [C II] in- 
tensity at this peak is I[cn] = (1-31 ± 0.15) x 10~ 4 
erg s _1 cm _2 sr _1 . This is in very good agreement with an 
earlier measurement of the integrated [C II] intensity by 
Stacey et al. (1991) who found J [CII] = (1.4 ± 0.3) x 10~' 4 



erg s _1 cm _2 sr _1 (Stacey 2001; see also Crawford et al. 
1985). A comparison of the FIFI [C II] results with our 
ISO [C II] data of M51 will be presented in a separate paper 
(Nikola et al. 2002). Preliminary reduction of our [C II] 
ISO measurement of M51 shows a lower intensity com- 
pared to the FIFI intensity, but roughly consistent with 
the ratio of the beam sizes. 

There is also a peak in the north-east (NE peak) of M51 
(at RA = 13 h 27 m 53.3 s , Dec = 47°29'03" (1950)) coincid- 
ing with an enhancement in the optical emission in the 
spiral arm. The integrated [C II] intensity at the NE peak 
is 7r CII ] = 9.2 x 10~ 5 erg s _1 cm _2 sr _1 or about 70% of 
the nuclear value. A third peak in the [C II] emission in 
the south-west (SW peak) of M51 (at RA = 13 h 27 m 38.5 s , 
Dec = 47°25'33" (1950)) coincides with a position of a 
spiral arm where enhanced optical emission can be seen. 
Here the integrated [C II] intensity is I[cu] = 5.6 x 10 -5 
erg s _1 cm~ 2 sr _1 . 

The NE and SW peaks are roughly symmetrically 
aligned with respect to the nucleus, and fall on local bright- 
ness peaks in the optical spiral arms. Other than this, 
there is only a weak indication for spiral structure in the 
[C II] map. There is an extension in the [C II] emission that 
follows the northern spiral arm past the NE [C II] peak, 
but the NE and SW peaks appear to be quite localized, 
essentially point sources to our 55" beam. 




Fig. 1. — Contour map of the [C II] line superimposed on an optical image of M51. The peak intensity is ^[cil] 



and the contour lines are in steps of 1.5 X 10 erg 
beam size (55" FWHM) of the [C II] observation. 



starting at 1.5 X 10 



1.3x10 4 erg s 'cm 2 sr 1 
Thc hashed circle represents the 



3.2. Comparison with FIR Continuum 

A comparison of the distribution of the [C II] emission 
with the distribution of the FIR emission shows very good 
agreement. Figure 2 shows the contour map of the inte- 
grated [C II] line intensity superimposed on a contour map 
of the 170 /im FIR continuum (Smith 1982) at a spatial 
resolution (49") similar to that of the [C II] map. The 
location of the main peak and the NE peak in both dis- 
tributions match each other very well. There is no peak 
towards the SW, but there is a ridge from the nuclear re- 
gion that touches our SW [C II] peak. It is likely that 
these are the same structures, but not coincident due to 
modest signal-to-noise ratios in both maps. 

The integrated [C II] luminosity of M51 is L[cii] = 
3 x 10 8 L Q or 1% of the FIR luminosity as measured by 
the IRAS satellite (3 x 10 10 L Q , Rice et al. 1988, scaled 
to 9.6 Mpc). This is at the high end of the £[cii]/£fir 
ratio observed in galactic nuclei(Crawford et al. 1985; 
Stacey et al. 1991), but consistent with the reported ra- 
tios for late type galaxies (NGC 6946, Madden et al. 1993; 
NGC 4038/4039, Nikola et al. 1998). We use the Smith 
(1982) map to calculate the i[cii]/^FiR ratio within our 
beam. To do this, we correct the Smith (1982) luminos- 
ity estimates (that include only the 80 /xm and 200 /im 
contiunuum) to include the shorter wavelength FIR con- 
tinuum as observed by IRAS. This correction amounts to 
a factor of 1.5 increase in the Smith (1982) Ltr map. The 
ratios range form 0.6% at the nucleus to 1.1% and 1.4% 



respectively at the NE and SW peaks (Tablel). 

3.3. Comparison with CO (1-^0) Observations 

As shown by Crawford et al. (1985) and Stacey et al. 
(1991) the combination of the integrated [C II] intensity 
with the integrated intensity of the 12 CO (1— >0) line is 
a powerful diagnostic tool for measuring the star forma- 
tion activity. To do this analysis we derived a CO (1— >0) 
contour map from the CO data of Lord & Young (1990) 
because their observations have been carried out with a 
spatial resolution (45") similar to our [C II] observation. 
However, their selection of positions does not cover the 
whole galaxy and we interpolated over the narrow gaps 
in their observation. Unfortunately their coverage shows 
a gap at the position of the NE peak. Figure 3 shows 
the superposition of the [C II] contour map on top of the 
derived contour map of the CO (1— >0) line. The peak in- 
tegrated intensity of the CO emission is Tco = 8.4 x 10~ 8 
erg s _1 cm _2 sr -1 . At the [C II] peaks we derived inte- 
grated intensities of the CO line of Iqo — 6.6 x 10~ 8 



erg s 'cm 2 sr l at the main [C II] peak, Ico = 1-8 x 10 8 
at the NE peak, and I co = 1.7 x 10~ 8 



erg s 'cm 2 sr l 



erg s cm sr at the SW peak. The peak of the CO 
emission is not at the nominal nucleus of M51 and does 
not coincide with the main peak of the [C II] emission. 
Another observation of the CO (1— >0) line by Scovillc & 
Young (1983) with the same spatial resolution but worse 
spatial coverage shows the peak of the CO emission in 
the center of the galaxy. The CO (1— >0) map of Garcia- 
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Fig. 2. — Contour map of the [C II] emission (white) superimposed on the contour map of the FIR continuum emission (black) at 170 fim 
(Smith 1982). The gray area outlines the observed field in the [C II] line. The contour lines of the FIR continuum are 9, 13, 18, 25, 35, 50, 
and 70% of the FIR peak (8 X 10~ 18 W m^Hz^sr- 1 ). The (0,0) position corresponds to RA = 13 h 27 m 46.4 s , Dec = 47°27'13" (1950). The 
white hashed circle represents the beam size (55" FWHM) of the [C II] observation. The spatial resolution of the FIR continuum map is 49". 



L[CU] [Lq] 
Lfir [Lq] 
L\cU]/Lf IR 



Table 1 
Observed [C II] and FIR luminosity 



Main [C II] Peak 
"TO 7 



3.2 x 

5.4 x 10 9 
0.6% 



NE Peak 
2.2 x 10 y 
1.9 x 10 9 

1.1% 



SW Peak 

1.4 x 10 Y 

1.0 x 10 9 

1.4% 



a ) Averaged over the whole galaxy M51 (CI « 2.85 X 10 6 sr). 

b ) From Smith (1982) (corrected to include flux < 80 /im using IRAS numbers) 



Total M51 a ' 

3.0 x 10 s 

3.0 x 10 10 

1.0% 



Burillo, Guelin, & Cernicharo (1993), with a spatial res- 
olution of 33", shows a local minimum at the center and 
the CO peak shifted by about (-20", -20") relative to the 
center of M51. This position is again in better agreement 
with the CO peak of Lord & Young (1990). Apart from the 
slightly different positions of the CO peak and the [C II] 
main peak both distributions are reasonable similar. The 
distribution of the CO emission is also elongated in north- 
east south-west direction and shows an extension toward 
the SW peak. The overall agreement suggests much of the 
observed [C II] emission arises in PDRs on the surface of 
molecular clouds. 

3.4. Comparison with Ha Emission 

Figure 4 shows the superposition of the [C II] contour 
map on the Ha contour map of van dcr Hulst et al. (1988). 
Despite the very different spatial resolution (8" for the Ha 
map) the maps are similar in the sense that the NE peak 
and the SW peak coincide with an accumulation of strong, 
discrete Ha knots in the spiral arms. There are three mod- 



erately strong Ha peaks located in the spiral arm in the 
south-east where we only detect weak [C II] emission. This 
might be the result of incomplete coverage of the full ve- 
locity dispersion at the edge of our array as described in 
section 2. Pixels showing this effect have not been taken 
into account for the data reduction. These "missing" pix- 
els created an incomplete image plane near the edge of our 
field of view and can therefore easily cause the wavy con- 
tour line. Apart from this, the superposition of the [C II] 
map on the Ha map shows nicely that the [C II] emission 
peaks coincide with strong H II regions and that the [C II] 
emission therefore traces mainly PDRs in star formation 
regions. However, weak [C II] emission is also seen in re- 
gions where there is almost no Ha emission, e.g. in the 
west-north-west region of M51. 

4. ORIGIN OF THE [C II] EMISSION IN M51 

4.1. Cold and warm neutral medium 

Comparison of the H I column density maps of Weli- 
achew & Gottesman (1973) and Tilanus & Allen (1991) 
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Fig. 3. — Contour map of the the [C II] emission (white) superimposed on a contour map of the 12 CO(l— >0) emission (black) derived from 
CO data from Lord & Young (1990). The contours of the CO line are 5, 10, 15, 20, 30, 40, 50 K km s _1 (corrected for beam efficiency). The 



(0,0) position corresponds to RA 
[C II] observation. 



13 h 27 m 46.4 s , Dec = 47°27'13" (1950). The hashed circle represents the beam size (55" FWHM) of the 



(Figure 5) with our [C II] map shows that the distributions 
do not agree. The H I column density has a minimum at 
the center of M51 and shows a ridge of higher column den- 
sity around the nucleus. On top of the H I column density 
ridge are several local column density maxima. The NE or 
SW peaks lie on this H I column density ridge but do not 
coincide with the local column density maxima. Although 
local H I column density maxima are in the vicinity of the 
NE and SW peak the overall morphology of the [C II] dis- 
tribution and the H I column denisty distributions do not 
agree. It therefore seems unlikely that much of the peak 
[C II] emission arises from diffuse neutral medium. 

Here we calculate the expected integrated [C II] intensity 
from the cold and warm atomic medium. For the calcu- 
lation we assume that C + behaves like a two-level system 
and that the [C II] radiation is optically thin. This gives: 
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(see also Madden et al. 1993, 1997) where h is the Planck 
constant, v is the frequency of the transition, A is the Ein- 
stein coefficient for spontaneous emission, 2.29 x 10~ 6 s -1 
(Nussbaumer & Storey 1981), and g u /gi — 2 is the ratio of 
the statistical weights in the upper and lower level. The 
first term in the equation above is due to excitation by 
hydrogen impacts while the second term is due to excita- 
tion by electron impacts. The critical density for collisions 



with atomic hydrogen, n cr .H, is deduced from the cooling 
function for collisional excitation of [C II] by atomic hy- 
drogen (Launay & Roueff 1977), and that for collisions 
with electrons, n C r,o ; is deduced from collision strengths 
for [C II] collisions with electrons (Blum & Pradhan 1992), 
and fitted as functions of temperature. We further assume 
a gas phase abundance of carbon relative to hydrogen of 
X c + = 3 x 1CP 4 , that all the carbon is ionized, a beam 
filling factor of $b = 1 , and that the same range of con- 
ditions for the cold neutral medium (CNM temperature 
T = 50 - 100 K, density n H = 50 - 200 cm" 3 and ion- 
ization fraction x c = 3 x 10 -4 ) and the warm neutral 
medium (WNM temperature T = 4000 - 8000 K, density 
riH = 0.5 — 5 cm -3 and ionization fraction x c = 3 x 10 -2 ) 
(Kulkarni & Hciles 1987; Kulkarni & Heiles 1988) of the 
Milky Way apply for M51 . We derived the H I column den- 
sity for our calculation from the gas surface density values 
of Lord & Young (1990) who used H I data of Rots et al. 
(1990) smoothed to a resolution of 45" for their estimate. 
The results are given in Table 2. 

Clearly, emission from the WNM is not important for 
most of the regions mapped in [C II] . Even the maximum 
values for the expected integrated [C II] intensity from the 
WNM are below the [C II] detection limit. 

A non-negligible contribution to the [C II] emission from 
the WNM is only possible in the north-west where ex- 
tended [C II] emission at the detection limit is visible. 
At this position the H I column density derived from 
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Fig. 4. — Contour map of the [C II] emission (white) superimposed on a contour map of the Ha emission (black) (van der Hulst ct al. 
The spatial resolution of the Ho map is 8" and the contours are 0.05, 0.25, 0.5, 0.75, 1.0, 1.5, 2.0, 3.0, 4.0 xl0~ 13 erg s _1 cm -2 . 



the data of Lord & Young (1990) reaches its maximum 
value of Nh « 2.5 x 10 21 cm~ 2 . The estimated contri- 
bution of the WNM in this region goes from a minimum 
of 2.2 x 10~ 6 erg s _1 cm _2 sr _1 or about 7% of the in- 
tegrated [C II] intensity (< 3x 10 -5 erg s _1 cm _2 sr _1 ) 
up to a maximum of 2.3 x 10~ 5 erg s~ 1 cm _2 sr~ 1 or 77% 
over the range of the physical conditions of the WNM. 
In this region, in the north-west, the contribution of the 
CNM to the integrated [C II] intensity can easily ac- 
count for the total observed intensity in the north-west 
region. Even the minimum contribution from the CNM 
(2.9 x 10~ 5 erg s^ 1 cm _2 sr _1 ) can account for all of the 
observed [C II] emission. At standard conditions (tih ~ 
100 cm" 3 , T ~ 80 K) the contribution from CNM to the 
[C II] emission is 4.6 x 10~ 5 erg s _1 cm _2 sr _1 which is more 
than 1.5 times the observed [C II] emission. 

At the [C II] peak positions the estimated integrated 
[C II] intensity arising from the CNM can also be rela- 
tively high. The estimated fraction of [C II] emission com- 
ing from the CNM at the NE and SW peak is between 



- 15% and - 36% (see Table 2). Although the CNM con- 
tribution is smaller at the position of the main [C II] peak 
it is still not negligible (8 - 12%; see Table 2). 

Since we assumed that all carbon in the CNM and WNM 
is in the form C + and that all of the observed H I column 
density is associated with the [C II] emission, the esti- 
mated values are upper limits. Therefore a contribution 
to the [C II] emission form the CNM is likely smaller than 
that called out above, so that at the nucleus , NE and SW 
peaks, the contribution may be small (<10 - 20%). How- 
ever, in the north-west region of M51 the CNM could be 
the dominant source of the [C II] emission. 



4.2. Ionized medium 

The estimate of the contribution of the ionized medium 
to the [C II] emission is very controversial. We will discuss 
four different approaches here. The easiest approaches are 
similar to the one we used to estimate the contribution of 
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Fig. 5. — Contour map of the [C II] emission (white) superimposed on a contour map of the H I column density (gray scale) (Tilanus & 
Allen 1991). The spatial resolution of the H I map is 20" and the contours are 3, 6, 10.5, 14.9 xlO 20 atoms cm -2 . Darker colors correspond 
to higher column densities. 



Table 2 
Observed and estimated integrated [C II] intensity from the CNM, the WNM, and the ELDWIM 



Main fC II] Peak 

=3 



'[Cii] in [erg s 'cm "sr 
NE Peak SW Peak 



Total M51 



a) 



Observed /[ 



[CII] 



N K [cm- 2 ] 



1.31 x 10" 



9.6 x W 



'[CII]min 



CNM 



b) 



-'[CII]std 

T b ) 

-<[CII]max 

J[CII]min 
f[CII]std 
-f[CII]max 

EM [pc cm" 6 ] 



1.0 x 10 -5 

1.6 x 10" 5 (12%) 

1.6 x 10- 5 (12%) 

8.4 x 10" 7 (0.6%) 

1.7 x 10" 6 (1%) 
8.7 x 10" 6 (7%) 

1387 



9.2 x 10"° 

1.3 x 10 21 

1.4 x 10- 8 (15%) 
2.2 x 10" 5 (24%) 

2.2 x 10" 5 (24%) 

1.1 x 10" 6 (1%) 

2.3 x 10" 6 (2.5%) 

1.2 x 10 -5 (13%) 

412 



5.6 x 10" 



1.2 x 10 2 



3.7 x 10" 



1.1 x 10 2 



1.3 x 10~ 5 (23%) 1.2 x 10" 5 (32%) 
2.0 x 10" 5 (36%) 1.8 x 10" 5 (49%) 
2.0 x 10 -5 (36%) 1.8 x 10 -5 (49%) 

9.7 x 10" 7 (3%) 

K (5%) 

(27%) 



1.1 x 10" 6 (2%) 
2.1 x 10" 6 (4%) 
1.1 x 10- 5 (20%) 

653 



2.0 x 10" 6 (5%) 
9.9 x 10" 6 '"- ' 



91 



ELDWIM ^ CII l min 

'[CH]p 



3.9 x 10" 5 (30%) 6.6 x 10" 6 (7.2%) 9.4 x 10" 6 (16.8%) 7.8 x 10" 6 (21%) 
(> 100%) 3.9 x 10" 5 (42%) (> 100%) (> 100%) 



a ' Averaged over the whole galaxy M51 (f2 « 2.85 x 10 6 sr). 

b ) Temperatures and densities used for CNM: 

min.: n H = 50 cm" 3 , T = 100 K; max.: n H = 200 cm" 3 , T = 50 K; standard: n H = 100 cm" 3 , T = 80 K; 

c ) Temperatures and densities used for WNM: 

min.: n H = 0.5 cm" 3 , T = 8000 K; max.: n H = 5 cm" 3 , T = 4000 K; standard: n H = 1 cm" 3 , T = 6000 K; 

d ) see text for minimum and maximum contribution 

e ) Fractions of observed integrated [C II] intensity are noted in parenthesis. 



the CNM and WNM. This leads to the following formula 

^cxp(-ai) 



I\ 



[CII] 



hvA 
47rn cri e 



_i+(i+^cx P (-ai))^ 
xI c+ $ B / S n H n ds 

where n cr . c « 49 cm -3 is the critical density of C + for col- 
lisions with electrons at 8000 K (Blum & Pradhan 1992), 

and L riB_n c ds is the emission measure (EM). 

Klein, Wielebinski, & Beck (1984) and van der Hulst et 
al. (1988) determined the flux density of the thermal ra- 
dio continuum from single-dish measurements with a beam 
size of 76" and from interferometric measurements with a 
spatial resolution of 8", respectively. The thermal flux 
density derived from the interferometric measurement is 
about 34% of that derived from the single-dish observation 
and originates mainly from prominent giant H II regions 
with sizes of about 500 pc and electron densities of about 
n c «l-2 cm" 3 . These H II clouds are distributed within 
the nuclear region of M51 and along the spiral arms. In 
contrast, the distribution of the thermal 2 cm radio con- 
tinuum derived from the single-dish observation does not 
show any spiral signature. 

Due to the big beam size of FIFI the main contributor 
of the ionized medium within the beam is most likely the 
extended low density warm ionized medium (ELDWIM). 
In our first approach we therefore estimate a possible [C II] 
emission from ionized gas using the single-dish measure- 
ments of the radio continuum of Klein, Wielebinski, & 
Beck (1984) assuming that all the thermal radio contin- 
uum arises from the ELDWIM. In this case the intensity 
of the thermal radio continuum is the same within the FIFI 
beam as within the single-dish beam. From the emission 
measure, derived from the thermal radio continuum, and 
assuming the low density limit (n c <C n„) and a beam fill- 
ing factor of unity we estimate the contribution of the EL- 
DWIM to the observed [C II] intensity from the above for- 
mula. The flux density of the thermal radio continuum at 
the NE and SW peak are taken from the 2 cm thermal con- 
tinuum contour map of Klein, Wielebinski, & Beck (1984; 
their Figure 7) . The intensity of the thermal radio contin- 
uum at 2 cm (T2th) at the NE peak is T 2 th ~ 2.5 mK. From 
this we get an emission measure of EM w 215 pc cm" 6 
(see Spitzer 1978 for the calculation) resulting in an ex- 
pected integrated [C II] intensity from ELDWIM at the 
NE peak of J [C ii] « 1.8 x 10" 5 erg s^on^sr -1 . This is 
about 20% of the observed integrated [C II] intensity. At 
the SW peak we get an intensity of the thermal radio con- 
tinuum at 2 cm of T 2 th ~ 4 mK from the contour map in 
Klein, Wielebinski, & Beck (1984). This leads to an emis- 
sion measure of EM w 344 pc cm" 6 and an integrated 
[C II] intensity of I[ai\ ~ 2.9 x 10~ 5 erg s _1 cm~ 2 sr _1 
which is about 52% of the observed integrated intensity. 
Within an area of 60" x 70" at the center of M51, Klein, 
Wielebinski, & Beck (1984) determined a thermal radio 
continuum at 2 cm of 11 mjy. Assuming again the same 
intensity within the FIFI beam as within this area we cal- 
culate an emission measure of EM « 1382 pc cm" 6 and an 
expected integrated [C II] intensity of I[cu] ~ 1-18 x 10" 4 
erg s _1 cm _2 sr _1 . This is about 90% of the observed in- 
tegrated [C II] intensity. The flux density of the ther- 
mal 2 cm radio continuum of the whole galaxy M51 is 
S 2 th « 55 mjy (Klein, Wielebinsky, & Beck 1984). This 



flux density was measured within an elliptical area with 
half axis of 5.75' and 4.93'. The corresponding emission 
measure is EM w 91 pc cm" 6 . This results in an ex- 
pected mean integrated [C II] intensity from H II regions 
of -^ [Cii] ~ 7-8 x 10" 6 erg s" 1 cm" 2 sr" 1 for the whole galaxy 
which is about 21% of the observed mean integrated [C II] 
intensity. 

In the nuclear region and the north-east region the in- 
terferometer map of the 6 cm radio continuum of van der 
Hulst ct al. (1988) reveals numerous strong sources. Thus 
the assumption that all the thermal radio continuum radi- 
ation in these regions arises from ELDWIM might be too 
simple. In our second approach we therefore assume that 
34% of the thermal radio continuum arises from clumpy 
H II regions as indicated in the radio continuum map of 
van der Hulst et al. (1988). Correcting for beam filling 
from the 76" beam to the 55" beam and assuming an elec- 
tron density of n c = 1.5 cm" 3 for the clumpy medium 
enhances the expected integrated [C II] intensity originat- 
ing from the clumpy medium and ELDWIM by a factor of 
1.26 compared to the first approach. 

Except for the NE peak and the total galaxy both ap- 
proaches result in large fractions of the [C II] emission 
attributed to ELDWIM reaching 100% for the nucleus. 
However, by assuming a beam filling factor of unity and 
the low density limit these values are upper limits. 

Another approach to estimate the [C II] emission from 
ELDWIM is to compare the integrated [C II] intensity with 
integrated [N II] intensity. The [N II] emission comes en- 
tirely from ionized medium and the critical densities of the 
[N II] emission lines are similar to the critical density of 
the [C II] . Treating the N + ion as a pure three level system 
and assuming an abundance ratio for X(C + )/X(N + ) it is 
possible to estimate the contribution of the ELDWIM to 
the [C II] emission from [N II] measurements. In the low 
density limit the expected ratio of the integrated [C II] in- 
tensity to the integrated [N II] 122/xm intensity from EL- 
DWIM is 



I\ 



[CII] 



121.898 A x 



:,s 



'[Nil], 122 



157.741 A 12 
2expH|) 

702 



,158 



A, 



c+ 



721+720 



A 



N+ 



121.898 2.29 x 10" 



310 

157.741 X 7.46 x 10" 6 X ~49~ 
1.98 3 x 10" 4 



2.3X10- 8 



;10- 8 +4.7xl0- 9 



1 x 10- 



«9 (1) 

where A158 and A\ 22 are the emission coefficients for 
spontaneous emission of the [C II] line and the 122 //m 
[N II] line, respectively (Nussbaumer & Storey 1981, Nuss- 
baumer & Rusca 1979), the 7y's are the collision rate coef- 
ficients derived from the effective collision strengths given 



by Lennon & Burke (1994), and n„ e and 



,158 



are the 



critical densities for collisions with electrons at a tem- 
perature of 10 4 K for the [C II] line and the 122 /an 
[N II] line, respectively. The contribution of the ELD- 
WIM to the observed [C II] emission is then given by 
(9 x / [ Niip22)/ / [°cii] rVed - For the comparison it is prefer- 
able to have the [N II] and [C II] intensities obtained with 



the same instrument. Using our ISO [N II] 122 /im and 
ISO [C II] data (Nikola et al. 2002) we obtain a contribu- 
tion of the ELDWIM to the observed ISO [C II] emission 
of (9 x 1.9 x 10~ 6 )/3.9 x 10~ 5 w 44% for the NE peak, 
(9 x 1.4 x 10- 5 )/7.2 x 10~ 5 w 175% for the nucleus, and 
(9 x 2.7 x 10~ 6 )/2.5 x 10~ 5 w 97% for the SW peak, where 
the intensities are in units of [erg s" 1 cm~ 2 sr -1 ]. These 
numbers are higher than the values derived above in our 
first two approaches. Note that the calibration of our ISO 
data is still uncertain. Therefore the individual ISO in- 
tensities should only be used with reservations. However, 
the ratio of the ISO intensities should not be affected by 
the calibration uncertainty. The relative statistical errors 
of the ISO [N II] intensities are « 30% at the NE- and 
SW-peak and w 5% at the nucleus, and the relative statis- 
tical errors of the ISO [C II] intensities are « 10% at the 
NE-pcak, s=s 4% at the SW-peak, and w 2% at the nucleus. 

From recent COBE FIRAS observation Petuchowski & 
Bennett (1993) estimated the morphology of the ionized 
medium in the Milky Way. They estimate that only a 
small fraction (17%) of the [N II] 122/xm emission actually 
arises from ELDWIM and that the majority (82%) of the 
[N II] emission comes from externally ionized cloud sur- 
faces. These externally ionized cloud surfaces have higher 
densities (n H = 150 cm~ 3 ) than the ELDWIM. If we as- 
sume the same conditions in M51 as in the Milky Way then 
the low density limit used in the above calculation might 
therefore not be appropriate. Assuming again a pure three 
level system for the N + ion and a pure two level system 
for the C + ion the expected intensity ratio of the 158 /im 
line to the 122 /im line in an ionized medium with a den- 
sity of 150 cm -3 and a temperature of 10 4 K decreases 
to I15S/I122 ~ 2.8. If we divide the [N II] emission into a 
part (17%) arising from ELDWIM and a part (82%) aris- 
ing from more dense externally ionized cloud surfaces and 
using our ISO data we obtain a much lower contribution of 
the ELDWIM to the observed [C II] emission: 7.2% for the 
NE peak, 30% for the nucleus, and 16.8% for the SW peak. 
For reference Petuchowski & Bennett (1993) estimate that 
the ELDWIM in the Milky Way contributes about 53% to 
the [C II] emission. This fraction is higher than the frac- 
tions we derive using the estimated [N II] morphology from 
COBE for all our positions but similar or smaller than the 
fraction we estimated for the SW peak and the nucleus in 
the first three approaches. For the NE peak this contribu- 
tion would be higher than the values derived above in all 
four approaches. 

For the minimum ELDWIM contribution to the [C II] 
emission we listed the estimated integrated [C II] intensity 
derived from the COBE FIRAS [N II] 122 /im morphol- 
ogy (Table 2). For the total galaxy M51 the only way to 
estimate the minimum value was by our first approach. 
All the maximum contribution from ELDWIM are higher 
than 100%, except for the NE peak. We got the maximum 
contribution to the [C II] emission at this position when 
we assumed that all the ISO [N II] emission arises from 
ELDWIM. 

In the further analysis we will use the minimum [C II] 
contribution from the ELDWIM. 



4.3. Photodissociation regions 



Strong [C II] emission is commonly attributed largely to 
PDRs associated with molecular cloud surfaces. However, 
as shown in the previous section in the case of M51 compo- 
nents other than PDRs can account for a large fraction of 
the observed integrated [C II] intensity. But presently it is 
not possible to make firm statements about how much of 
the [C II] emission is really coming from other components 
of the interstellar medium other than PDRs. Therefore we 
will follow four ways in analyzing the [C II] emission from 
PDRs. First we assume that all the observed integrated 
[C II] intensity originates from PDRs ("case-1"). Then 
we subtract the minimum contribution of all non-PDR 
contributions from the observed integrated [C II] inten- 
sity ("case-2"). As a third case we subtract the standard 
contributions of the CNM and the WNM and the minium 
contribution of the ELDWIM from the observed [C II] 
intensity and attribute the remaining intensity to PDRs 
("case-3"). For the last case we subtract the maximum 
contribution from CNM, WNM, and ELDWIM from the 
observed integrated intensity ("case-4"). In this case the 
NE peak is the only position where a PDR contribution 
to the [C II] emission is left over. 

4.3.1. Determination of density, FUV intensity, and 
beam filling factor of PDRs 

To determine the density, the far-UV (FUV) intensity 
and the beam filling factor of the PDRs we compare the 
integrated [C II] intensity with the integrated 12 CO(1^0) 
intensity from Lord & Young (1990) and the intensity of 
the FIR continuum from Smith (1982) as corrected us- 
ing IRAS data and relate this to PDR model predictions 
(Wolfire ct al. 1989, Stacey et al. 1991). We determine 
these values for the three points of interest, the main [C II] 
peak, the NE peak, and the SW peak. In addition we de- 
termine the mean values for the entire galaxy M51. The 
values used in this estimate are given in Table 3. 

Figure 6 shows the densities (solid lines) and the 
FUV intensities (dashed lines) as functions of Y[cu] = 
J[Cii]/Xfir and Y~ co = /co/xfir from PDR models 
(Wolfire, Hollenbach, & Tielens 1989, Stacey et al. 1991), 
where xfir is the integrated FIR intensity expressed in 
units of the local Galactic interstellar radiation field \o = 
2 x 10~ 4 erg s _1 cm _2 sr _1 (Draine 1978). Using the ratios 
Yjcn] and Too eliminate the contribution of the beam fill- 
ing factor if it assumed that the [C II] and 12 CO (1 — ► 0) 
line and FIR continuum emission are associated with the 
same PDRs. The solid lines connect PDR models with the 
same density while dashed lines connect PDR models with 
the same incident FUV intensity. At a density of about 
10 4 cm -3 the Vjcii] ratio reaches a maximum. This results 
in a degeneracy of the PDR model solutions. If the (lco> 
yjcii]) point falls above the 10 2 cm~ 3 density curve then 
two different solutions (low and high density) exist for the 
density and the FUV intensity. The filled symbols repre- 
sent "case-1" (all [C II] from PDRs), the star symbol rep- 
resent "case-2" (minimum contribution of [C II] emission 
of non-PDRs subtracted) and the open symbols represent 
"case-3" (standard contribution of [C II] emission of non- 
PDRs subtracted) for the individual positions. The gray 
bar shows the range of possible [C II] contributions origi- 
nating in PDRs. If all the [C II] emission would arise from 
non-PDRs then the gray bar reaches I'jcii] = 0. The (Yco> 
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Table 3 
Paramaters used for the PDR model 





Main [C II] peak 


NE peak SW peak Total M51 a) 


/co b,c) 

XFIR d) 


6.6 x 10"* 
112 


1.8 x 10" 8 1.7 x 10"* 2.8 x 10"* 
40 20 23 


T b) 
-'[CII] 

-f[cii]/-fco 


1.31 x 10" 4 
2000 


Observed 
9.2 x 10" 5 5.6 x 10" 5 3.7 x 10" 5 
5000 3300 1300 


T b) 
1 [CII] 


Minimum Contribution from non-PDRs subtracted 
8.0 x 10" 5 7.0 x 10" 5 3.1 x 10" 5 1.6 x 10" 5 


r b) 

-'[CII] 


Standard Contribution from the CNM and the WNM 

and minimum of ELDWIM subtracted 
7.3 x 10" 5 6.0 x 10" 5 2.3 x 10" 5 8.5 x 10" 6 


r b) 

-'[CII] 


Maximum Contribution from non-PDRs subtracted 
1.7 x 10" 5 



a ) Averaged over the whole galaxy M51 (CI Ri 2.85 x 10 6 sr). 

' The integrated intensities are in units of [erg s cm" sr" 1 ]. 
c ) 7 co = 1.6 X 10" 9 / (ip/T/c) di/ [erg s^on^sr- 1 ] (from Lord & Young 1990) 
d ' From Smith (1982) (corrected to include flux < 80 ^m using IRAS numbers); 
in units of the local interstellar radiation field: 



\o = 2 X 10 4 erg s 



- 1 (Draine 1978) 



yjcii]) positions of "case-1" for the NE peak and the SW 
peak are just at the edge of the range for a possible PDR 
solution but fall clearly within the range of a solution if 
the uncertainty of the [C II] measurement is taken into 
account. 

The PDR solutions for the main [C II] peak, the NE and 
SW peak, and the total galaxy that we derived from this 
plot arc shown in Table 4. As can be seen from the PDR 
plot the FUV intensities for all the positions in M51 are 
modest. The beam filling factor which can be determined 
by the intensity ratio of the observed FIR radiation and 
the modeled FUV radiation (Wolfire, Hollcnbach, & Tie- 
lens 1989, Stacey et al. 1991) is also shown in Table 4. The 
calculation of the beam filling factor assumes that about 
all the FUV radiation impinging on the PDR is converted 
into FIR continuum. We use this ratio instead of the ratio 
of the observed to the predicted [C II] intensity because of 
the uncertainty of the fraction of [C II] originating from 
PDRs. Also the change of the modeled FUV intensities is 
small for the different [C II] contributions from PDRs. For 
the area at the main [C II] peak we get a very high beam 
filling factor of nearly 100% for the low density PDR solu- 
tion. In all other cases the beam filling factor is between 
7% and 38%. 



4.3.2. Mass and column density of the [C II] emitting 

regions 

We estimate the column density and the mass of the 
[C II] emitting photodissociated gas using the results of 
the PDR model above. For the calculation we assume op- 
tically thin emission. The column density and mass is then 



given by: 



V H 



47rA/i 



[CII] 



hcAX c+ & b 



-expl 



-f) 



1+f exp(-ai)- 



a 



cxp(- 



and Mi 



H 









M 



O 



where Q is the area in steradian, M p is the proton mass, 
Mq is the mass of the sun, D is the distance, h is 
the Planck constant, c is the vacuum speed of light, 
A is the Einstein coefficient for spontaneous emission, 
Xq+ = 3 x 10" 4 is the abundance of ionized carbon, 
g n and g\ are the statistical weights for the upper and 
lower level, x c « 3 x 10" 4 is the ionization fraction, and 
ftcr.H = 3.5 x 10 3 is the mean critical density for collisions 
with atomic and molecular hydrogen. The density and the 
beam filling factors were taken from the PDR results (Ta- 
ble 4). In this calculation we assumed that all carbon is in 
singly ionized form and that the gas temperature is 200 K. 
For this gas temperature we assume a critical density for 
collisions with electrons of n C r,c ~ 10 cm" 3 . This critical 
density is deduced from the collision strengths given by 
Blum & Pradhan (1992). 

The results for the column density and the mass are pre- 
sented in Table 5. For comparison, the molecular masses 
at NE and SW peak and at the nucleus within a FIFI beam 
and of the total galaxy M51 are also determined. The mass 
of the molecular gas was derived from the column density 
of the molecular hydrogen given in Lord & Young (1990). 
The mass of the [C II] emitting gas averaged over M51 and 
at the nucleus is between 2 - 34% of the molecular mass. 
At the NE and the SW peak the fraction of the mass of 
the [C II] emitting photodissociated gas is as high as 72%. 

5. DISCUSSION 
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Fig. 6. — PDR-plot, showing the density — FUV-fiux phase space of PDRs as a function of ^jcill = ^[Clll/XFIR and Yqq = ^co/xfir- 
The solid lines and the dashed lines show PDR models at fixed densities and fixed FUV intensities, respectively. Peak [C II] observations 
and an average over M51 arc represented by a filled symbol. The location of the filled symbol in the plot marks the phase space for the 
case that the total observed integrated [C II] intensities originate from PDRs. The gray bars show the range in phase space that can be 
covered by varying the contribution of [C II] emission originating from PDRs for each observed position. The stars represent the case that the 
minium contribution from non-PDRs is subtracted form the observed integrated [C II] intensity. The open symbols represent the case with 
the standard contributions from non-PDRs subtracted from the observed [C II] emission and the cross represents the case with the maximum 
contribution from non-PDRs subtracted from the observed integrated [C II] intensity. 



Table 4 
Results from the PDR model 





Main 


C II] peak 


NE peak SW peak 


Total M51 a) 






Density 


Density Density 


] 


Density 






Low 


High 


Low High Low High 
From observed values 


Low 


High 


riH [cm 


"I 


300 


1.2 x 10 5 


w 10 4 900 3 x 10 3 


200 


6.3 x 10 5 


b ) 
Xfuv 


[Xol 


130 


460 


300 130 200 


60 


200 


$B C) 




0.86 


0.24 


0.13 0.15 0.10 


0.38 


0.12 








Minimum non-PDRs contribution subtracted 




riH [cm 


"I 


^ 100 


4.5 x 10 5 


750 2 x 10 4 270 1 x 10 5 


— 


5.4 x 10 5 


Xfuv 


[Xol 


100 


520 


350 420 100 300 


— 


250 


$B C) 




1.12 


0.22 


0.11 0.10 0.20 0.07 


— 


0.09 








Standard non-PDRs contribution subtracted 




riH [cm 


- 3 l 


^ 100 


4.6 x 10 5 


700 3 x 10 4 140 1.7 x 10 5 


— 


1.4 x 10 6 


Xfuv 


[Xol 


100 


530 


200 500 80 300 


— 


250 


$B C) 




1.12 


0.21 


0.20 0.08 0.25 0.07 


— 


0.09 








Maximum non-PDRs contribution subtracted 




riH [cm 


- 3 l 


— 


— 


1 x 10 6 


— 


— 


Xfuv 


[Xo] 


— 


— 


760 


— 


— 


$B C) 




— 


— 


0.05 


— 


— 



a ) Averaged over the whole galaxy M51 (H pb 2.85 X 10 

b ) In units of the local interstellar radiation field: 
Xo = 2 X 10~ 4 erg s _1 cm _2 sr _1 (Drainc 1978) 



sr). 



c ) Beam filling factor: <3? B = XFIR 
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Table 5 
Column density and mass of the [C II] emitting gas compared to the molecular mass 



Total M51 a ) 



Main [C II] Peak 



NE Peak 



SW Peak 



M v 



[M ] 



7V H [cm- 2 ] 
M H [M ] 

N m [cm" 2 ] 
M m [Mo] 

N m [cm- 2 ] 
M m [M ] 



AT H i [cm- 2 ] 
M H i [Mq] 



1.0 x 10 9 

Density 
Low High 



2.4 X 10 21 
1.5 x 10 8 



7.6 x 10 2() 
4.6 x 10 7 



3.1 x 10 
1.9 x 10 s 



-' 4.8 x 10 20 
2.9 x 10 7 



2.8 x 10 21 
1.7 x 10 8 



4.6 x 10 2() 
2.8 x 10 7 



2.6 x 10 8 

Density 
Low High 



2.5 X 10 8 

Density 
Low High 



From Observed Values 
1.5 x 10 21 2.5 x 10 21 
8.9 x 10 7 1.5 x 10 s 



1.8 x 10 21 
1.1 x 10 s 



Minimum non-PDR contribution subtracted 
4.8 x 10 21 1.2 x 10 21 2.7 x 10 21 6.2 x 10 20 
2.9 x 10 8 7.4 x 10 7 1.6 x 10 s 3.8 x 10 7 

Standard non-PDR contribution subtracted 
2.4 x 10 21 1.2 x 10 21 2.9 x 10 21 4.5 x 10 20 
1.4 x 10 8 7.4 x 10 7 1.8 x 10 s 2.7 x 10 7 

Maximum non-PDR contribution subtracted 
4.5 X 10 20 
2.7 X 10 7 



1.35 x 10 1() 

Density 
Low High 



2.2 X 10 21 4.1 X 10 20 
4.6 x 10 9 8.5 x 10 8 



2.3 x 10 2() 
4.9 x 10 8 



1.2 x 10^ J 
2.6 x 10 8 



a ) Averaged over the whole galaxy M51 (CI RJ 2.85 X 10 6 sr). 

b ' The molecular mass was derived from the H2 column densities given in Lord & Young (1990). 



5.1. Starbursts at the NE and SW Peaks 

The [C II]/ 12 CO(1^0) line intensity ratio is an indica- 
tor of OB star formation activity in dusty galaxies (Stacey 
et al. 1991). A ratio of 4400 is typically found in starburst 
nuclei, while a value near 1200 is more typical of quiesent 
spiral galaxies. In the north-east of M51, near the NE 
peak, we obtain the highest [C II] to CO integrated inten- 
sity ratios in M51 (5000; for "case-1") similar to starburst 
nuclei. From their data Lord & Young (1990) also derived 
the highest star formation efficiency in the north-east re- 
gion close to the NE peak. The moderate integrated in- 
tensity ratio at the nucleus (I[cn]/Ico ~ 2000) indicates 
moderate star formation activity. Although the peak of 
the [C II] emission is located at the center of M51 there is 
also a large amount of CO in this region resulting in the 
moderate intensity ratio. In addition the nucleus of M51 
is most likely an AGN (Ford et al. 1985; Makishima et al. 
1990; Kohno et al. 1996; Grillmair et al. 1997) and not 
a starburst region. The integrated intensity ratio in the 
south-west of M51 (I[cu]/Ico ~ 3300) is smaller than at 
the NE peak. Although we can not give a firm value for 
the fraction of [C II] emission from PDRs in M51 at the 
NE peak at least 50% or more of the integrated [C II] in- 
tensity originates most likely from PDRs. This result and 
the high [C II]/CO(l— >0) intensity ratio suggest that the 
most active star forming region in M51 lies in the north- 
east region, probably near the NE peak. The maximum 
of the [C II] to CO integrated intensity ratio, however, is 
next to the NE peak between the spiral arms. This is likely 
the result of interpolation in the CO map. No CO data 
were obtained by Lord & Young (1990) at the exact loca- 
tion of the NE peak. The low [C II]/CO(1-»0) intensity 
ratios at the nucleus and the SW peak may be the effect 
of beam dilution. If the ISM in these regions consists of 
only a few but highly active star forming knots within an 
extended, less active region the [C II]/CO(l— >0) intensity 
ratio averaged over the FIFI beam would be low. Such a 



morphology of the ISM is suggested e.g. for NGC 4038/39 
(Nikola et al. 1998). 

The finding of enhanced star formation activity in the 
north-cast and probably also in the south-west is also 
supported by the distribution of the Ha emission (van 
dcr Hufst et al. 1988), the FUV emission (Bcrsicr et 
al. 1994), and the 12 CO(2^1) emission (Garcia-Burillo, 
Guclin, & Cernicharo 1993), the FIR continuum (Smith 
1982), and the thermal radio continuum at A = 2 cm 
(Klein, Wiclcbinski, & Beck 1984). All the distributions 
show emission which either peaks in the north-east and 
south-west or only peaks in the north-east but is elongated 
in north-east south-west direction. However, this distribu- 
tion is most pronounced in the [C II] emission with a main 
focus towards the north-east region. 

Additional support for the enhanced star formation ac- 
tivity and the special location of the activity comes from 
analysis of the kinematics and dynamics of M51. Tully 
(1974 c) concluded in his analysis that the density-wave 
pattern in M51 is confined between the inner Lindblad 
resonance and corotation and that the spiral arms outside 
corotation can be accounted for by material clumping in- 
duced by the galaxy interaction. Corotation occurs at a 
radius of about 5.8 kpc (Tully 1974 c; converted to a dis- 
tance of M51 of 9.6 Mpc). Both the NE peak and the 
SW peak lie at about the same distance from the center 
of M51 as corotation. Thus they are located at the transi- 
tion zone of the density-wave pattern and the spiral arms 
which are attributed to material clumping. Especially the 
SW peak coincides with a position of extreme clumping 
(Tully 1974 c; position "n" in Figure 10 therein). This is 
also the position where the orbits of all the test particles 
which were originally at or beyond 6.7 kpc pass through 
(Toomre & Toomre 1972). Material clumping can also be 
seen at the position of the NE peak. The encounter model 
of Toomre & Toomre (1972) shows crowding of test par- 
ticles along the outer spiral arm which starts at the NE 
peak where the density-wave arm terminates and follows 
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on to north and west to the southern tidal tail. This spiral 
arm can be attributed to material clumping (Tully 1974 
c). The relatively weak [C II] emission that we observe in 
the north-west region might originate from this clumped 
material. The coincidence of the location of the NE and 
SW peak with the region where most or all trajectories of 
the test particles pass through and where the density wave 
arms terminate might indicate that the enhanced [C II] 
emission is caused by an enhanced mass flow crossing the 
spiral arms. This could result in a high rate of cloud-cloud 
collisions which then trigger star formation. The position 
of the NE and SW peak at these special locations where 
the tidal tials emerge from the disk in the simulations of 
Toomre & Toomre (1972) also indicates that the star for- 
mation there is triggered by the galaxy interaction. The 
trajectories of test particles crossing the spiral arms shown 
by Toomre & Toomre (1972) are almost concentrated in a 
single spot in the south-west in contrast to the north-east 
where the crossing area is more spread out. This would 
result in a beam dilution effect at the SW peak, resulting 
in a smaller [C II]/CO(l— >0) intensity ratio. 

5.2. Comparison with other Spiral Galaxies 

Except for the NE and SW peak and the main peak in 
the center of M51 the [C II] emission is distributed very 
uniformly over the whole galaxy. There is no clear signa- 
ture of the spiral structure visible in the [C II] emission. 
This could be due to a smoothing effect given our beam 
size or the presence of an underlying [C II] emission in 
M51 or both. However, [C II] observations with FIFI of 
NGC 6946 (Madden et al. 1993), which is at about the 
same distance as M51 but a little bit less extended, shows 
some spiral structure. Also the [C II] map of M83 taken 
with FIFI (Gcis ct al. in preparation), which has a similar 
extension as M51 reveals very clumpy [C II] emission fol- 
lowing the spiral arms. This difference in the NGC 6949 
and M83 maps and the M51 map supports the presence of 
an underlying extended [C II] emission in M51. We have 
shown in the previous sections that the CNM and the EL- 
DWIM can in principle contribute a large fraction to the 
[C II] emission. The ELDWIM is most likely the main con- 
tributor to an underlying [C II] emission in most of M51. 
But in the north-west region of M51 the CNM is most 
likely the main source of the [C II] emission. How much 
the ELDWIM contributes to the [C II] emission in M51 
remains unclear. Even with the use of the [N II] 122/im 
line it was not possible to present a firm statement about 
the real contribution of the [C II] emission from ionized 
gas. Mapping of M51 in the [C II] and both [N II] lines 
with high spatial resolution might help to disentangle the 
contributions from the various ionized gas phases and to 
distinguish between an on-arm, inter-arm, and underly- 
ing [C II] emission. High spatial resolution would also be 
required to find out if the PDRs lie on a separate spiral 
arm (maybe between the ionized arm and the molecular 
arm or coinciding with the molecular arm) and to distigu- 
ish between [C II] emission from a possible PDR arm and 
emission from neutral gas in the H I arm. 

The PDRs at the NE and SW peak might be located in 
molecular superclouds of masses -Mmsc ~ 10 7 M Q within 
the giant molecular association with masses of Mqma ~ 
10 7 — 10 8 M Q observed by Vogel, Kulkarni, & Scovillc 



(1988) and Rand & Kulkarni (1990) in these regions. This 
is supported by the low FUV intensity of xfuv ~ 10 2 Xo 
which we estimated from the PDR models. Such low FUV 
intensities arc normally expected in PDRs in giant molec- 
ular clouds (M GM c « 10 5 - 10 6 M Q ) (Stacey ct al. 1991, 
1993) and might therefore also be the case in PDRs in 
molecular superclouds. 

In the nucleus the high density PDR solution is sup- 
ported by the detection of the HCN molecule in the center 
of M51 and in the inner spiral arms (Nguyen et al. 1992; 
Kuno et al. 1995; Kohno et al. 1996). This finding sug- 
gests that the density of the molecular gas in these clouds 
is > 10 5 cm~ 3 . No observations of this density tracer 
are available in the outer parts of M51 and therefore no 
statement can be made about which PDR solution is more 
likely 

6. SUMMARY 

We obtained a map of M51 in the [C II] fine structure 
line at 158 /xm with a spatial resolution of 55" (FWHM). 
The [C II] emission arises from the whole visible extent 
of M51. The main peak of the [C II] emission is located 
at the center of M51 and has an integrated intensity of 



I\ 



[CII] 



1.31 x 10 4 erg s 'cm 2 sr 1 . 



In addition we 



found a [C II] peak in the north-east and in the south-west 
of M51. The NE peak is stronger than the SW peak and 
has an integrated [C II] intensity of 70% of the integrated 
intensity of the main peak. The NE peak and the SW peak 
lie roughly on a straight line going through the center of 
M51 and they have the same distance to the center as coro- 
tation of the density wave pattern. Both peaks are located 
at the position where the density wave pattern terminates 
and where the spiral structure is attributed to material 
clumping due to the interaction of M51 with NGC 5195. 
The [C II] emission in the NE peak originates mainly from 
PDRs. We also find that the strongest star formation ac- 
tivity occurs in the north-east of M51 close to the NE peak. 
We suspect that this enhanced star formation activity is 
triggered by cloud-cloud collisions due to a high mass flow 
crossing the spiral arms. There might also be high star 
formation activity in the south-west triggered by the same 
effect but which is hidden to our observation due to beam 
dilution. 

From the PDR model we derive two possible solutions 
for the individual positions we investigated in M51. The 
estimated FUV intensity is similar for both solutions with 
an intensity of a few 100 times the local Galactic interstel- 
lar radiation field. For the density we obtain very differ- 
ent solutions with n ~ 10 2 — 10 3 cm" 3 for the low density 
solution and n ~ 10 4 — 10 6 cm~ 3 for the high density so- 
lution. At the center of M51 the high density solution is 
supported by the detection of HCN in this region. The 
PDRs in the NE and SW peak are most likely associated 
with the molecular superclouds and giant molecular asso- 
ciations found in these regions. This is supported by the 
relatively low FUV intensity derived from the PDR model. 
The total mass of the [C II] emitting photodissociated re- 
gions in M51 is M « 2.6 x 10 8 M©. This is only a small 
fraction (2%) of the total molecular mass. If all [C II] 
arises in PDRs that fraction can be as high as 34%. In the 
nucleus the fraction of the [C II] emitting photodissociated 
gas ranges between 3% and 19% and at the NE and SW 
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peak the mass of the [C II] emitting gas can be as high as 
72%. 

A large fraction of the overall [C II] emission in M51 can 
originate in an underlying extended medium. It remains 
unclear how much the ELDWIM contributes to this emis- 
sion, but it could be the main source. The CNM can also 
be a significant source of [C II] emission. Especially at the 
north-west part of M51 where most of the [C II] emission 
might originate from CNM. At the NE and SW peak a 



smaller fraction of the [C II] emission might emerge from 
CNM. 
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